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Abstract 



Our goal in this paper is to examine the discovery potential of laboratory ex- 
periments searching for the oscillation v^fe) — > u T , in the light of recent data 
on solar and atmospheric neutrino experiments, which we analyse together 
with the most restrictive results from laboratory experiments on neutrino os- 
cillations. In order to explain simultaneously all present results we use a 
four-neutrino framework, with an additional sterile neutrino. Our predictions 
are rather pessimistic for the upcoming experiments NOMAD and CHORUS, 
which, we find, are able to explore only a small area of the oscillation param- 
eter space. On the other hand, the discovery potential of future experiments 
is much larger. We consider three examples. E803, which is approved to 
operate in the future Fermilab main injector beam line, MINOS, a proposed 
long-baseline experiment also using the Fermilab beam, and NAUSICAA, an 
improved detector which improves by an order of magnitude the performance 
of CHORUS/NOMAD and can be operated either at CERN or at Fermilab 
beams. We find that those experiments can cover a very substantial frac- 
tion of the oscillation parameter space, having thus a very good chance of 
discovering both v, L — > v T and v e — ► v T oscillation modes. 
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I. INTRODUCTION 



Present data from solar and atmospheric neutrino experiments favour the hypothesis of 
neutrino oscillations. Nevertheless, this interpretation requires confirmation from further 
experiments, in particular from laboratory experiments. All solar neutrino experiments |l ] 
find less v e than predicted theoretically. However, the uncertainties in the calculation of the 
solar neutrino flux may still be rather large || and new experiments ||,f|] are being planned 
to further investigate the possible origin of the solar neutrino deficit. As for atmospheric 
neutrino experiments, two of them P-[71 measure a ratio v^/v e smaller than expected from 
theoretical calculations. Here, in addition to the uncertainties on the estimation of the 
atmospheric neutrino flux and on the neutrino- nucleon cross section |§, one has to add the 
uncertainties due to the modest data sample. New experiments are also being planned in 
this area ||]. 

A complementary approach is to look for neutrino oscillations in laboratory experiments, 
where the experimental conditions, in particular the shape, energy, and flux of the neutrino 
beam are under control. A number of neutrino experiments have recently started taking data. 
These are the LSND experiment at Los Alamos ||, which looks for z/ M — > v e oscillations, and 



the CHORUS [JTOj] and NOMAD [[□]] experiments at CERN, which are primarily sensitive to 
v ji{ye) — > v t oscillations. Recent data from LSND may be consistent with the existence of 



neutrino oscillations [[C|, although no formal claim has been made so far by the collaboration. 
First results from CHORUS and NOMAD should be available in 1996. In addition, a 
number of new v^iv^ — > v T oscillation experiments are being discussed. At Fermilab, 
a new, very intense v^{v e ) beam is planned to be available around the year 2001. Two 
complementary experiments are being discussed to benefit from this beam. These are E803 



||13|| , a short-baseline experiment similar in design to CHORUS, and MINOS ||14j| , a long- 
baseline experiment, which proposes to detect the neutrinos produced at Fermilab with 
a 10 Kton detector located in the Soudan mine, around 730 km away from the neutrino 
source. Furthermore, new experimental techniques are being studied to upgrade NOMAD 
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and CHORUS [15,|16|. Ultimately, these new techniques could result in the design of new 



detectors able to improve the performance of NOMAD and CHORUS/E803 by one order of 



magnitude P7,I5|,I5|. Such detector(s) could be operated either at CERN or at Fermilab 
beams. 

In this paper we examine the prospects of success of all these experiments in the light 
of present data. CHORUS and NOMAD were conceived at a time when the dominant 
scenario for neutrino masses was consistent with very light is e , and a v T of about 10 
eV constituting the hot component of the dark matter. This scenario arose as the most 
natural solution to explain the solar neutrino problem, simultaneously providing a candidate 
for the hot dark matter. However, present data no longer favour this simple scenario. 
Indeed, explaining the results from solar and atmospheric neutrino experiments in the usual 
three-neutrino framework is very difficult and requires a large degree of fine-tuning. One 
has to choose between throwing away part of the data and considering a larger scheme. 
The "minimal" scheme to explain all data without fine-tuning seems to be a four-neutrino 
framework (u e , u T , v s ) where v s is a sterile neutrino. Using this framework we perform 
a consistent analysis of the data from solar and atmospheric neutrino experiments as well 
as results on neutrino-oscillation laboratory experiments. This analysis enables us to re- 
assess the discovery potential of CHORUS and NOMAD as well as to study the prospects 
for the new experiments discussed above. We find that CHORUS and NOMAD have a 
rather marginal chance of discovering v^{v e ) — > v T oscillations. However, the future v T 
oscillation experiments have much better prospects since they cover a very large fraction of 
the oscillation parameter space. 

The outline of the paper is as follows. In Sec. || we review the formalism for neutrino 
oscillation in a general multi-family framework. Section [TIT] is devoted to a summary of 
the present experimental status for solar and atmospheric neutrino experiments as well as 
accelerator and reactor neutrino oscillation experiments. The basic ingredients of the four- 
neutrino framework are presented in Sec. [TV] and the results of our analysis of the available 
experimental data, discussed in Sec. |T|, are displayed in Sec. |V} Section |V| describes 



succinctly the upcoming and future v T oscillation experiments while Sec. [VI]| is dedicated 
to the study of the prospects for the discovery of v T oscillations. Finally we present our 
conclusions in Sec. IV 111. 



II. FORMALISM 

If neutrinos have a mass, the weak eigenstates v a produced in a weak interaction (i.e., 
an inverse beta reaction or a weak decay) will be, in general, a linear combination of the 
mass eigenstates i/j 

n 

v a = J2U ai Vi (1) 
1=1 

where n is the number of neutrino species and U is a unitary mixing matrix. 

After travelling a distance L, the neutrino can be detected in the charged-current (CC) 
interaction v N' — > lp N with a probability 

n n /A- A 

P a p = ^-4EE MUaiU^Upj] sin 2 ( ) (2) 
i=i j=i+i \ £ J 

The probability, therefore, oscillates with oscillation lengths A^- given by 

A^ = |m?-mj| L/g = Am?- L/E 
2 ' eV 2 m/MeV ' eV 2 m/MeV 1 J 

where E is the neutrino energy. Each experimental set up has a different characteristic value 
of the ratio L/E and is thus most sensitive to oscillation lengths such that Am^ ~ 1/(L/E). 
The typical values of L/E are summarized in Table |. 

In general neutrino beams are not monochromatic. Thus, rather than measuring P a p, 
the experiments are sensitive to the average probability 



(Pap) 



J dE u —a cc {E v )P a(3 {E v )e{Ev) 

J dE ^ a cc(EMEu) (4) 

n n /A- A 

^-4E E MUaiU^U^isin 2 ^ ), 

i=l j=i+l \ ^ / 
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where $ is the neutrino energy spectrum, acc is the cross section for the process in which 
the neutrino is detected (in general a CC interaction), and e(Ej) is the detection efficiency 
for the experiment. For oscillation lengths such that Am| ^> 1/(L/E) the oscillating phase 
will have been over many cycles before the detection and therefore it will have averaged to 
(sin 2 (Aj,-/2)) = 1/2. On the other hand, for Am| <C 1/(L/E), the oscillation did not have 
time to give any effect and (sin 2 (Aj :) /2)) « 

For the case of two-neutrino oscillations the above expressions take the well known form 

Pap = S al 3 - (26 a(S - 1) sin 2 (20) (sin 2 (^\) (5) 



since 



U 



( cos 9 sin 9 ^ 
y — sin 9 cos 9 i 



(6) 



Most neutrino oscillation experiments present their result in the two-family mixing lan- 
guage as regions in the plane (Am 2 , sin 2 (26 1 )). Using the previous expression it is possible 
to translate their results into transition probabilities. 



III. EXPERIMENTAL STATUS 

A. Accelerator and Reactor Neutrino Experiments 

There are two types of laboratory experiments to search for neutrino oscillations. In 
a disappearance experiment one looks for the attenuation of a neutrino beam primarily 
composed of a single flavour due to the mixing with other flavours. At present the most 
restrictive experiments of this kind are the reactor experiment at Bugey fll9 |, which looks 



for v e disappearance, and the CDHSW experiment |20| at CERN, which searches for 
disappearance. Both experiments show no indication of neutrino oscillation. Their results 
are presented as exclusion areas in the two-neutrino oscillation approximation in Fig. |I[ We 
translate their results into limits on the transition probabilities: 
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(P ee ) > 0.93 from Bugey for Am 2 > 4eV 2 

(7) 

(P w ) > 0.95 from CDHSW 

Both results are given at 90% CL. The maximum probability from Bugey is larger for 
smaller mass differences since the neutrino flux normalization can be better determined by 
the experiment. However, in the range of masses we are interested the relevant limit is the 
one given above. 

In an appearance experiment one searches for interactions by neutrinos of a flavour not 
present in the neutrino beam. The most restrictive experiments are the E776 experiment 



at BNL |21] for the 9^ — > u e channel and the E531 experiment at Fermilab [^] for the 
— > v T channel. Neither of these experiments shows evidence for neutrino oscillation 
on those channels. Their results are also presented as exclusion areas in the two-neutrino 
oscillation approximation in Fig. |I[ In terms of transition probabilities we find 

(P eM ) % 1.5 x 10- 3 from E776 

(8) 

(P^r) < 2 x 10~ 3 from E531 

at 90% CL. 

Recently the Liquid Scintillator Neutrino Detector (LSND) experiment ]I2] has an- 
nounced the observation of an anomaly that can be interpreted as neutrino oscillations in 
the channel — > i> e . Most of the oscillation parameters required as explanation are already 
excluded by the E776 [^TJ and KARMEN experiments. Their result can be marginally 
compatible with these previous results for Am 2 w 6-8 eV 2 and mixing sin 2 (2#) ~ 3 x 10~ 3 . 
One must point out, however, that these ranges are still quite far from certain. 



B. Solar Neutrinos 

At the moment, evidence for a solar neutrino deficit comes from four experiments [|TJ. 
Putting all these results together seems to indicate that the solution to the problem is not 
astrophysical but must concern neutrino properties. The standard explanation for this deficit 
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would be the oscillation of v e to another neutrino species either active or sterile. Different 
analyses have been performed to find the allowed mass differences and mixing angles in 
the two-flavour approximation ||24|| . They all seem to agree that there are three possible 



solutions for the parameters: 

• vacuum oscillations with Am^ = (0.5-1) x 1CT 10 eV 2 and sin 2 (26>) = 0.8-1 

• non-adiabatic-matter-enhanced oscillations via the MSW mechanism [^5| with Am^ = 
(0.3-1.2) x 10~ 5 eV 2 and sin 2 (2#) = (0.4-1.2) x 10" 2 , and 

• large mixing via the MSW mechanism with Am 2 ei = (0.3-3) x 10~ 5 eV 2 and sin 2 (26>) = 0.6- 
1. 

C. Atmospheric Neutrinos 

Atmospheric neutrinos are produced when cosmic rays (primarily protons) hit the at- 
mosphere and initiate atmospheric cascades. The mesons present in the cascade, decay 
leading to a flux of v e and which reach the Earth and interact in the different neutrino 
detectors. Naively the expected ratio of to v e is in the proportion 2:1, since the main 
reaction is 7r — >■ fiv^ followed by /i — > ev^v e . However, the expected ratio of muon-like inter- 
actions to electron-like interactions in each experiment depends on the detector thresholds 
and efficiencies as well as on the expected neutrino fluxes. 

Currently four experiments have observed atmospheric neutrino interactions. Two ex- 
periments use water-Cherenkov detectors, Kamiokande |5]J| and 1MB j7|, and have observed 
a ratio of z^-induced events to z/ e -induced events smaller than the expected one. In par- 
ticular Kamiokande has performed two different analyses for both sub-GeV neutrinos || 
and multi-GeV neutrinos [HI, which show the same deficit. On the other hand, the results 



from the two iron calorimeter experiments, Frejus |26] and NUSEX ||27|| , appear to be in 
agreement with the predictions. 

The results of the three most precise experiments are shown next. They are given as 



a double ratio R^/ e /R^ of experimental-to-expected ratio of muon-like to electron-like 



S 



events. The expected ratio Rff/^ is obtained by Monte Carlo calculation of the atmospheric 
neutrino fluxes. We have used the expected fluxes from Gaisser et al. [^j (see also Ref. 
p9| ) and Volkova |3(J depending on the neutrino energies (see discussion in Sec. |TV| ). Use 
of other flux calculations plj would yield similar numbers. 



Rn /e /R%° = 0.55 ±0.11 for 1MB 

R v/e/Rf?/e = °- 6 ± °- 09 for Kamiokande sub-GeV ^ 
R p/e/R%e = °- 59 ± 0- 1 for Kamiokande multi-GeV 
R v/e/Rf?/e = 1.06 ± 0.23 for Frejus 
The most plausible explanation for this anomaly is to suppose that oscillates into another 
flavour. The oscillation — > v e is almost completely ruled out by the reactor experiment 



data [JL{J. We are then left with — > v T oscillations or oscillations to a sterile neutrino. 
The allowed range of masses and mixings in the two-family approximation from a global fit 
to the previous data is shown in Fig. 0.a and can be summarized as 

Amjj > 2 x 10~ 3 sin 2 (2#) = 0.56-1 . (10) 

The experiments have also observed an angular dependence on the value of these double 
ratios. The previous ranges were obtained without using the angular information. An 
analysis based on the angular dependence of the ratio leads to an upper limit on the possible 
value of Am 2 , as can be seen from the Kamiokande multi-GeV analysis and the 1MB 
data 0. However the allowed regions obtained from the best fit to these two sets of data 
are inconsistent at the 2a level (see Fig. 5 in Refs. [§,0). For this reason, we have chosen 
not to use the constraints arising from the angular dependence of the data in this analysis. 

D. Dark Matter 

There is increasing evidence that more than 90% of the mass in the Universe is dark and 
non-baryonic. Neutrinos, if massive, constitute a source for dark matter. Stable neutrinos 
can fill the Universe of hot dark matter if their masses add up to a maximum of about 30 



eV. However, scenarios with only hot dark matter run into trouble in the explanation of 
the formation of structures on small scales of the Universe. On the other hand, models for 
structure formation favour the presence of cold dark matter. These models, however, fail in 
reproducing the data on the anisotropy of the microwave background. Currently, the best 
scenario to explain the data considers a mixture of 70% cold plus 30% hot dark matter . 
This translates into an upper limit on neutrino masses [j33|| : 

<4r-7eV. (11) 

i 

IV. FOUR-FLAVOUR MODELS 

Naive two-family counting shows that it is very difficult to fit all experimental information 
in Sec. [TTT] with three neutrino flavours pl| , p6| , even without invoking the LSND data. The 



solar neutrino deficit could be due to u e — > oscillations and the atmospheric neutrino 
deficit to —>■ v T oscillations with the appropriate mass differences, for example with a 
mass hierarchy m„ r ^> m v ,m Ue . However, fitting this together with the present laboratory 
limits leaves no room for hot dark matter since the maximum allowed mass is of about 



vh Vt < 0.7 eV |29[. The only possible way out is to require that all three neutrinos are 
almost degenerate. This requires a certain degree of fine-tuning in order to explain the 
neutrinoless double beta decay data |36| , |37| . [] 



One could have z/„ — > v T oscillations for the atmospheric neutrino deficit with almost 



degenerate and v T with masses m v = m Vr fa 2.5 eV and m Ve fa ||36| , |39[| , but leaving 
out the explanation for the solar neutrino deficit. Or m v = m Vr fa eV and m Ue fa 2.5 to 
explain the atmospheric data but leaving unexplained both solar neutrino deficit and dark 



matter pSfl . 



1 Notice that this scenario is also inconsistent with the oscillation parameters observed by LSND, 
should this experiment confirm their results 
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Also, it is possible to explain the solar neutrino deficit with u e — > v T (u) with almost 
degenerate v e and v T (p) with masses m„ e = 77V, , ~ 2.5 eV and tti u , t) ~ 0, but leaving the 
atmospheric neutrino deficit unexplained p7| , |38[] . Also m„ e = Trip.. ~ eV and fn v , T) ~ 2.5 



eV would explain the LSND data if confirmed but leaves both atmospheric and dark matter 
without explanation. 

In the spirit of Pauli, one is tempted to introduce a new neutrino as a "desperate solution" 
40fl to understand all present data. The nature of such a particle is constrained by LEP 



results on the invisible Z width as well as data on the primordial 4 He abundance. Those rule 
out the existence of additional, light, active neutrinos. In consequence the fourth neutrino 
state must be sterile. 

There are different mass patterns that one can construct with four such neutrinos to 
verify all the experimental constraints and evidence for neutrino oscillation and masses. We 
will assume a natural mass hierarchy with two light neutrinos with their main projection in 
the v s and v e directions and two heavy neutrinos with their largest component along the 
and v T flavours. We will also require that the sterile neutrino does not mix directly with the 
two heavy states. As we will see, this is necessary to verify the constraints from big bang 



nucleosynthesis PTJ. Such a hierarchy appears naturally, for instance if one advocates an 
L e ±L M =i=L r discrete symmetry for the mass matrix [|34}| - In Ref. |36| a similar mass pattern is 
also generated via a combination of see-saw mechanism and loop mechanism. Other possible 



patterns which explain all data such as some inverted mass schemes |37| require a somehow 
large degree of fine-tuning to explain the absence of neutrinoless double beta decay. 

For any matrix presenting this natural mass hierarchy, the mixing matrix can be 
parametrized in a general way as 



2 A similar mass matrix is also introduced in Ref. 
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"5 jiT^efj,CeT C^iT^er 







(12) 



with Cj = cos^j and Sj = sin0j. For the sake of simplicity we have assumed no CP violation 
in the lepton sector. In this approximation m 1; m 2 <C rn 3 , m 4 and \m\ — m\\ <C m§, m\. Also, 
us and z/4 constitute 30% of dark matter in the Universe. This implies 7713 ~ 777.4 = 2-3.5 eV. 
Such a mass pattern has been argued to yield satisfactory results in Cold+Hot Dark Matter 



scenarios [331. We can define 



Am\ T 



AM 2 DM 



m; 



rrin 



12 21 
\rri3 — 77I4I 



(13) 



(4 - 10) eV 2 



Transition probabilities between the different flavours will now have contributions from 
the three oscillation lengths due to the three different mass differences in the problem which 
we will denote sin 2 (A so ; ar /2), sin 2 (A J 4 T /2), and sin 2 (A DM /2), respectively. We are now 
ready to reanalyse the experimental data presented in Sec. [Hi] in the four-flavour framework 
considering the oscillations with the three different oscillation lengths. 



V. GLOBAL ANALYSIS 

A. Laboratory Experiments 

Laboratory experiments are insensitive to oscillations due to the solar mass difference 
Am 2 solar . Also, for most of them the oscillation due to the dark matter mass difference will 
be averaged to 1/2. As a consequence the negative results presented in Sec. |TT^ will impose 
severe constraints on the mixing angles. 
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For the Bugey reactor experiment the relevant transition probability is the v e survival 
probability. For any value of the atmospheric mass difference this probability will always 
verify 

0.93 < < 1 - 24,4.(1 - 4/4) =► c\/ er > 0.96 . (14) 

For CDHSW the relevant probability is the survival probability 

0.95<P^ DH5W <l-^sin 2 (2^) sin 2 (2^ eM ) < 0.1 . (15) 

For E776 the situation is somehow more involved, since the value of the oscillating phase 
(sin 2 (A£i A f /2)) varies in the range A DM = 4-10 eV 2 due to the wiggles of the resolution 
function of the experiment (see Fig. [l]). Also, the experiment is sensitive to the atmospheric 
mass difference: 

1.5 x 10- 3 > P e f 76 = sin 2 (2# eM )c 2 T sin 2 (^p) 

+ g sin(2V) sin(20 eM ) sin(2# er ) cos(2^ T )c e/i - sin 2 (2^ r )( S 2 M c 2 r - s 2 r )c 2 M ] sin 2 (^) . 

(16) 

For any value of the atmospheric mass difference and the /xr mixing angle, the previous limit 
is verified if 

sin 2 (2^)c 2 r < (2-5) x 10~ 3 (17) 

The limit from E531 on the mixings e/x and er is always less restrictive than the previous 
ones for any value of Am 2 AT and AM|, M . 

Combining these constraints we obtain that e/x and er mixings are constrained to 

sin 2 (2^) < (2-5) x 10- 3 ^ 
sin 2 (2# er ) < 0.16 

where the range of sin 2 (26^) depends on the specific value of AMp M . 

If we now turn to the effect due to the oscillation with A at, we can rewrite the relevant 
probabilities for the different experiments expanding in the small angles e/x and er: 
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PCdhsw „ x _ i sin 2( 2 ^) - sin 2 (2^ T ) sin 2 (^) 

P E531 „ sin 2( 2 ^ r ) c 2 T si n 2 (^) (19) 
P e 7 ; 6 ~ sin 2 (2^)c 2 T sin 2 (^) + sm 2 (26 llT )sl sin 2 (^) . 

With the constraints in Eq. (Pf), the Bugey experiment is not sensitive to oscillations with 

In Figs. @ and [3] we show the exclusion contours in (Am 2 AT , sin 2 (2^ T )) which are due to 
the three experiments for different allowed values of the e/i and er mixings. 



B. Solar neutrinos 



Different authors have considered the propagation of neutrinos in the framework of more- 



than-two-neutrino oscillations [29,42-44|. Rather than reanalysing the solar neutrino data 



we will follow and adapt the results in Ref. [29]. Following the approach of Ref. [43 



we 



will express the transition probabilities in the framework of four neutrinos in terms of the 
two-family ones. 

For the solar neutrino deficit the relevant transition probability is the v e survival proba- 
bility. For solar neutrino distance and energies, both oscillations with the Am AT and AM|> M 
will have averaged to 1/2 and the survival probability in vacuum takes the form 

psoiar „ i _ 2c 2 M c 2 r (l - c 2 ^) - c\/ er sin 2 (2^) <sin 2 (^)> , (20) 

where we have expanded in the small mixings e/i and er. Given the constraints on the 
product c 2 M c 2 r (Eq. (0)), the effect of the new mixings is very small and the vacuum 
oscillation solution remains basically unchanged. 



For propagation in matter we follow the notation given in Ref. [fj3| . Following a similar 
procedure we can write the transition probability in the four-neutrino scenario in terms of 
the transition probability for two neutrinos as 

-^ee = 1 — C efi C er + C e/i C er-^ee J (21) 
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where is computed in the two-family case from the evolution equation in matter substi- 



tuting the electron density N e by N e c t \j? er - According to the results of Ref. [p^j , and taking 



into account Eq. (0), the two MSW solutions are still valid in the presence of the new mix- 
ings. The effect of the mixings appears to go in the direction of favouring the small mixing 
solution |T2U] . The large mixing solution is also in conflict with the constraints from big bang 



nucleosynthesis. As can be seen from Fig. 1 and Fig. 2 in Ref. for v e -v s oscillations this 
solution would lead to an excess on the primordial 4 He abundance excluded by the present 
data. 

C. Atmospheric Neutrinos 

There are in the literature several analyses of the atmospheric neutrino data in terms of 
neutrino oscillations in two-family and three-family scenarios [^,^,^], the last ones being 



mostly in the one-mass- dominance approximation. This approximation is not valid in our 
scheme since the atmospheric neutrino fluxes can show the effect of oscillations due to two 
oscillation lengths A at and Adm- We reanalyse the data in Eq. (pi), taking into account 



the effect of these two oscillations. In doing so we will follow the notation of Ref. ||46|| . 

In each experiment the number of /i events, N^, and of e events, N e , in the presence of 
oscillations will be 



where 



N„ = <tP m ) + Nl(P v ) , JV e = «P 6e > + N a IP^) , (22) 



^ = / lE^lk <E * )dWiC0S ^ <23) 



and 



<^> = 4 / ii^ p ^ik eiE " )dEMos s -> ■ (24) 

Here E v is the neutrino energy and is the flux of atmospheric neutrinos u a ; Ep is the 
final charged lepton energy and e{Ep) is the detection efficiency for such charged lepton; a 
is the interaction cross section v N — > N' I. The expected rate with no oscillation would be 
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N° 

pA/C _ (cyr\ 



The value of this ratio is different for each experiment as it depends on the threshold for the 
detected lepton energy as well as on the detection efficiency for the different lepton flavours. 
The double ratio R^/ e /R^P of the expected ratio of muon-like to electron-like events with 
oscillation to the expected ratio without oscillations is given by 



We perform a global fit to the data in Eq. (|9|) using the neutrino fluxes from Ref. [28| for 



E u < 3 GeV. For E v > 10 GeV we use the fluxes from Ref. [!J[J, and for energies in between 
the two fluxes are smoothly connected. We perform the integrals starting at the different 
thresholds of each experiment and include the published detection efficiencies |)H7|,[26H . For 



the neutrino interaction cross section the dominant process is the quasielastic cross section 
47fl . For larger neutrino energies (for Frejus and Kamiokande multi-GeV analyses) one-pion 



and multipion cross sections become relevant and are included j|8 . 

The results of our x 2 fit are shown in Figs. |2| and [3]. In Fig. |2] the results are shown 
for zero mixings e/i and er as in a two-family scenario. Figure |3] shows the effect of the 
inclusion of the mixings. As seen in the figure the inclusion of the er mixing leads to a more 
constrained area for the oscillation parameters. When mixing e\x and er are non-zero and 
taking into account the constraints in Eq. flT8p, Eq. (EH) takes the form 



^ ~ 1 + \ sin 2 (2M - sin 2 (2V)(sin 2 (^)> , (27) 

where we have used the fact that for atmospheric neutrino experiments (sin 2 (A so ; ar /2)) = 
and (sin 2 (A DM /2)) = 1/2. For the purpose of illustration we have used the approximation 
iVl/ 'N® = N® e /N® ~ 0.5. The effect of the er mixing in Eq. ( |27D is to increase the value 
of the double ratio. Therefore a larger amount of /ir oscillation is needed to account for the 
deficit. Due to the small values allowed, a non-zero mixing e/x does not modify the analysis 
of the atmospheric neutrino data provided that (AM|, M , sin 2 (26 , eM )) are allowed by the E776 
data (see Fig. [I]). 
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VI. u T OSCILLATION EXPERIMENTS 



A. CHORUS and NOMAD 

The two upcoming v^{v e ) — > v T experiments, CHORUS and NOMAD, are v T appearance 
experiments, i.e., they search for the appearance of v T 's in the CERN SPS beam consisting 
primarily of z/^'s, with about 1% u^s, as shown in Fig. [| (solid lines) |]10|,|lTJ. The mean 
energy of the beam is around 30 GeV and the detectors are located approximately 800 
m away from the beam source. The v T contamination of the SPS beam is virtually zero. 

CHORUS and NOMAD have complementary techniques to identify a r lepton. The NO- 
MAD experiment |TT| distinguishes a v T CC interaction from ordinary or v e interactions 



by exploiting the fact that the r lepton produced in a v T CC interaction, decays emitting one 
or more neutrinos which result in a measurable amount of missing transverse momentum, 
in the general direction of the charged lepton. NOMAD is in essence an electronic bubble 
chamber, with a continuous target alternating panels of a light material and drift chambers, 
followed by a transition radiation detector and an electromagnetic calorimeter. All the above 
detectors are located inside a 0.4 T magnetic field. Thus NOMAD measures essentially all 
the charged tracks and photons in the event, which enables a good reconstruction of the 
transverse missing momentum magnitude and direction. CHORUS ]10[] on the other hand, 



seeks to observe the finite path of the r on its emulsion target. At SPS energies, the r 
mean decay length is about 1 mm, giving two distinctive signatures in CHORUS emulsion. 
One-prong decays result in a short track with a kink, while three-prong decays appear as a 
short track splitting in three. Thus CHORUS rely on purely vertex techniques to identify 
the t. However, in order to reduce the number of events to be scanned, loose kinematical 
cuts are also applied. To this extent, the emulsion target is followed by a spectrometer and a 
compensating calorimeter. In the surviving events, tracks reconstructed in the spectrometer 
are extrapolated to the emulsions in order to determine where to scan. 

The initial goal of the SPS is to deliver 2.4 x 10 19 protons on target over two years to both 
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experiments. The intense neutrino beam allows very light detectors (NOMAD has a fiducial 
mass of near 3 t, while the emulsion target of CHORUS is 800 kg). The main characteristics 
of NOMAD and CHORUS as well as their expected performance are summarized in Table 

E 

B. Future u T Experiments at CERN 

There are a number of future z/ M (z/ e ) — > v T experiments beyond CHORUS and NOMAD 
being discussed at present. In the more immediate future there exists the possibility of an 
extended run at the CERN SPS beam, spanning several years after the initial period. Several 
suggestions have been made to upgrade CHORUS and NOMAD. The upgrade proposed 
for CHORUS would go in the direction of substituting the emulsion target with an active 



target made of scintillating capillaries read-out by CCDS [01. NOMAD, on the other hand, 



could be upgraded by adding an instrumented target made of a sandwich of passive, low-Z 



material (carbon or beryllium) and silicon detectors Ultimately, one would like a precise 
measurement of the event vertex, which, when combined with the event kinematics, results in 
a much improved sensitivity. Specifically, the sensitivity of the upgraded NOMAD detector 



could be improved by one order of magnitude [Tj|. A different experimental approach, 
suggested in pT7f] , proposes a 100 t liquid-CH 4 TPC in a high magnetic field (2 t). The idea 
is to use the quasi-free protons provided by the CH4 to completely constrain the kinematics 
of the z/ T -proton collisions. One can determine the momentum of the outgoing neutrino 
in the case of a r decay involving one neutrino, and reconstruct invariant masses. The 
t signal would appear as a mass peak over a very low background reduced by cuts a la 
CHORUS/NOMAD. 

As a specific example of this improved detector we have considered a suggestion to 
upgrade the NOMAD detector (see |T5| where detailed calculations are carried out). We 
will refer to this future detector with the generic name of Neutrino Apparatus with Improved 
CApAbilities (NAUSICAA). Such a device would permit r detection using both kinematical 
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and vertex techniques. Neutrino interactions would occur in a target made of a sandwich of 
silicon detector and a light- Z material. The event vertex would be precisely measured using 
the silicon information. The target would be followed by several planes of drift chambers 
and a calorimeter, with the full detector inside a magnetic field of around 1 T. The detector 
performance is summarized in Table [TV[ NAUSICAA can effectively suppress the CC and 
Neutral Currents (NC) v^Ve) backgrounds, thanks to the combination of two independent 
signatures of the r, i.e., its "long" lifetime and its decays with neutrino emission. 

Finally, there are several proposals for long-baseline experiments using the CERN SPS 
beam (see for example |^9| and references therein). 



C. Future v T Experiments at FNAL 

At Fermilab, a neutrino beam will be available when the main injector becomes opera- 
tional, around the year 2001. Compared with the CERN SPS beam, the main injector will 
deliver a beam 50 times more intense, but with an average energy around one third of that 
of the SPS neutrinos |]|[[J (see Fig. |). 

There are currently two experiments proposed to operate in this beam. One is a short- 
baseline experiment, E803 |13|], very similar in design and fiducial mass to CHORUS. Nev- 
ertheless, E803 is foreseen to have a sensitivity ten times better than CHORUS, thanks 
essentially to a much larger data sample (in addition to a more intense beam, the experi- 
ment expects a 4- year run). MINOS Jl4| is a long-baseline experiment, which proposes two 



detectors, separated by 732 km. This experiment can perform several tests to look for a 
possible oscillation z/ M — > v T in the small mass difference range. Of those, the most sensi- 
tive is the comparison of the fraction of CC-like (defined for v^) and NC-like events (a v T 
event, appears, except for the decay r — > [iv^Vr as a neutral current in the Minos detector). 
The main characteristics and expected performance of E803 and MINOS are summarized in 
Table |TTJ. 

Finally we will consider the possibility of installing the NAUSICAA detector (see Table 
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IVD as an alternative or a successor to E803 in the Fermilab beam. 



VII. DISCOVERY POTENTIAL 

We now turn to the study of the prospects for the discovery of v T oscillations in the 
upcoming CERN experiments CHORUS and NOMAD, in the future experiment E803 at 
Fermilab as well as a possible improved experiment, NAUSICAA, operated both at CERN 
and Fermilab beam. As an example of long-baseline experiment we also analyse the prospects 
for MINOS. 

After implementing the limits derived in Sec. [V] and considering the sensitivity of the 
experiments, one sees that for all facilities the only observable v e — > v T transition oscillates 
with an oscillation length A^m such that 

P er ^sin 2 (2Msin 2 (^) (28) 

Figure [| shows the regions accessible to the experiments in the (sin 2 (2# eT ), AM|, M ) plane. As 
can be seen, the almost complete accessible range of parameters for the CERN experiments 
CHORUS and NOMAD is already excluded by the Bugey data in the favoured range for 
the mass difference due to dark matter considerations. E803 at Fermilab and NAUSICAA 
at CERN could observe these oscillations in the whole favoured range of AM|, M , provided 
that sin 2 (26 l eT ) > 3 x 10 -2 . Furthermore, NAUSICAA operated at the Fermilab beam can 
go as low as sin 2 (2# er ) > 4 x 10~ 3 . MINOS, being a disappearance experiment is not 
sensitive to this oscillation. 

For transitions — > v T a four-neutrino framework predicts (unlike the naive two-family 
framework) two oscillations, dominated by the characteristic lengths A^m and A at. All 
experiments are in principle sensitive to both oscillations depending on the values of the 
mixing angles 

P° M ~ sin 2 (2£ eM ) sin 2 (£ er ) sin 2 (^) 

(29) 

Pj? ~ sin 2 (2^ T ) cos 2 (# er ) sin 2 (^) 
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In Fig. we show the regions accessible in the (sin 2 (2# eT ), AM^ M ) plane to the different 
experiments for an optimum value of the er mixing sin (29 eT ) = 0.16. As seen in the figure, 
the whole parameter space accessible to the CERN experiments CHORUS and NOMAD is 
already ruled out by the E776 data. E803 and NAUSICAA at CERN can marginally see 
this oscillation in the favoured range for AMp M for maximum er mixing. In particular if 
the LSND data are confirmed, this oscillation could be observable at both experiments for 
this optimum value of the er mixing. For smaller values of the er mixing these oscillations 
became very marginal or invisible. NAUSICAA at Fermilab can observe this oscillation for 
sin 2 (2# eAt ) sin 2 (# er ) > (0.3-1) x 10 -4 . If the LSND results are confirmed, the same oscillation 
could then be detected in NAUSICAA at Fermilab, provided sin 2 (6' er ) > 10~ 2 . As seen in 
the figure the accessible range for MINOS is already excluded by the E776 data. 

Figure [7| shows the region accessible to the experiments in the (sin 2 (20 MT ), Am AT ) pa- 
rameter space for different values of the other mixings. As seen in the figure, the possibility 
of observing this oscillation depends on the values of the er and e/i mixings. For zero er 
mixing, the oscillation can be observed by all experiments in the region allowed by all labora- 
tory and atmospheric neutrino experiments (see Fig. |7].a). It must be pointed out, however, 
that CERN experiments can observe this oscillation very marginally only. For larger values 
of the er mixing the situation is almost the same as long as the e/i mixing remains zero. If 
both mixings take their maximum allowed value the whole parameter space for the CERN 
experiments is ruled out by the atmospheric and E776 experiments data (see Fig. [7].b). 

E803 and NAUSICAA at CERN can observe oscillations with Am\ T > 5 x 10~ 2 eV 2 for 
any allowed value of the er and e/i mixings. NAUSICAA operated at the Fermilab beam 
could access oscillations with Am 2 AT > (2x 10~ 2 ) eV 2 . []. MINOS can reach mass differences 



3 There is even the possibility |5C] that NAUSICAA can be operated at two points, near the 
beam and at a few tenths of km allowing a limit Am 2 AT > (3 x 10 -3 ) eV 2 sufficient to close the 
atmospheric window 
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as low as Am 2 AT > (10 3 ) eV 2 , closing down the allowed window for the atmospheric neutrino 
experiments. 

VIII. CONCLUSION 

In this paper we have studied the discovery potential of laboratory experiments search- 
ing for the oscillation f M (z/ e ) — ¥ v T , after considering existing data on solar and atmospheric 
neutrino experiments as well as the results from laboratory experiments on neutrino os- 
cillations. To understand all these results in a common framework, we have introduced a 
fourth sterile neutrino and we have performed a comprehensive reanalysis of the data in 
this four-neutrino scenario. As an outcome of this analysis we can predict the number of 
expected events at future v T oscillations experiments at CERN and Fermilab for the allowed 
oscillation parameters. We find that at these facilities u e — > v T oscillations are dominated by 
mass differences of the order of AM|, M ~ 4-10 eV 2 while — > v T oscillations could occur 
with both AM 2 DM and Am 2 AT ~ 0.3-10" 3 eV 2 . 

Our predictions are rather pessimistic for the upcoming experiments NOMAD and CHO- 
RUS, which, we find, are able to explore only a small area of the oscillation parameter space 
in both u e — > u T and — > u T channels. The prospects are much better for future v T ex- 
periments. E803 and/or an improved detector (NAUSICAA) at CERN can explore v e — > v T 
with sin 2 (2# er ) > 10~ 2 . As for — > u r , the oscillations with AM|, M could be observed 
at E803/NAUSICAA-CERN in the range favoured by the LSND results, while oscillations 
with Am AT are accessible for mixing values allowed by the atmospheric neutrino data if 
Am 2 AT > 5 x 10- 2 eV 2 . 

The proposed long-baseline experiment MINOS at Fermilab has the unique potential of 
completely exploring the region dominated by Am AT thus, if present atmospheric neutrino 
data are correct MINOS must observe the oscillation — > v T . Finally, a future high- 
performance detector such as NAUSICAA at Fermilab could have the potential of exploring 
a very sizeable region of the parameter space for all the oscillation modes. The u e — > u T , 
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oscillations could be detected in the limit sin 2 (26> er ) ~ 1CT 3 and AM|, M ~ 10 eV 2 ; — > v T 
oscillations dominated by AM|, M could be detected in the limit sin 2 (2# eM ) ~ 3 x 10~ 4 and 
AMp M ~ 10 eV 2 , and — > z/ r oscillations dominated by Am\ T could be detected in the 
limit Am 2 AT ~ 2 x 10~ 2 (~ 3 x 10 -3 ) and maximal mixing. Indeed, NAUSICAA has the 
potential to observe both v„ — > v T and u e — > z/ r oscillations. 



ACKNOWLEDGMENTS 

We want to thank P. Lipari for providing us with the atmospheric neutrino fluxes. We 
acknowledge discussions with R. Vazquez, P. Hernandez, R. Sundrum and S. R. Mishra. J. 
J. G. -C. acknowledges the hospitality of Harvard University. 



23 



REFERENCES 



[1] GALLEX Collaboration, Phys. Lett. B327, 377 (1994); SAGE Collaboration, Phys. 
Lett. B328, 234 (1994); Homestake Collaboration, Nucl. Phys. B38 (Proc. SuppL), 47 
(1995); Kamiokande Collaboration, Nucl. Phys. B38 (Proc. SuppL), 55 (1995). 



[2] 
[3] 

K 

[5 

[6] 
[7] 



[9 

[io; 
[ii 

[12 
[13 
[14 



J. N. Bahcall and M. H. Pinsonneault, Rev. Mod. Phys. 64, 885 (1992). 

Borexino Collaboration, Gran Sasso National Laboratory Report LNSG-Report 94/99 
Vol I (1994); Icarus Collaboration, Nucl. Phys. B28A (Proc. SuppL), 486 (1992). 

Y. Totsuka, SUPERKAMIOKANDE, ICCR-Report 227-90-20 (1990); G. T. Ewan, Sud- 
bury Neutrino Observatory Proposal SNO-87-12 (1987). 

Kamiokande Collaboration, H. S. Hirata et al., Phys. Lett, bf B205, 416 (1988) and 
Phys. Lett. B280, 146 (1992). 

Kamiokande Collaboration, Y. Fukuda et al, Phys. Lett. B335, 237 (1994). 

1MB Collaboration, D. Casper et al, Phys. Rev. Lett. 66, 2561 (1991); R. Becker-Szendy 
et al, Phys. Rev. D46, 3720 (1992). 



P. Lipari, M. Lusignoli, and F. Sartogo, Roma preprint 1072-1994 and [hep-ph/9411341 



See for example W. C. Louis representing the LSND Collaboration, Nucl. Phys. B38 
(Proc. SuppL), 229 (1995). 

CHORUS Collaboration, N. Armenise et al., CERN-SPSC/90-42 (1990). 

NOMAD Collaboration, P. Astier et al., CERN-SPSLC/91-21 (1991), CERN- 
SPSLC/91-48 (1991), SPSLC/P261 Add. 1 (1991). 

H. Whyte, "Neutrino oscillations and neutrino mass", talk, MIT, 6 April 1995. 
P-803 Proposal, Fermilab, October 1993. 
P-875 Proposal, Fermilab, March, 1995. 

24 



[15] J. J. Gomez-Cadenas, A. Bueno and J. A. Hernando, A Proposal to Upgrade the NO- 
MAD Detector NOMAD/Memo 95-013(1995) 

[16] S. Buontempo et aJ., CERN-PPE/94-142 (1994), submitted to Nucl. Instrum. Methods; 
G. Martellotti et al, Nucl. Instrum. Methods A315, 177 (1992). 

[17] C. Rubbia, CERN-PPE/ 93-08 (1993), Proceedings of Conversaciones de Madrid, 5 
Sept., 1992, El Escorial, Madrid. 

[18] K. Winter, CERN-EP/89-182, Proceedings of the Protvino UNK Workshop (1989), 
p.83. 

[19] B. Achkar et al., Nucl. Phys. B424, 503 (1995). 

[20] CDHSW Collaboration, F. Didak et al., Phys. Lett. B134, 281 (1984). 

[21] E776 Collaboration, L. Borodvsky et al., Phys. Rev. Lett. 68, 274 (1992). 

[22] E531 Collaboration, Phys. Rev. Lett. 57, 2898 (1986). 

[23] B. Armbruster et al, Nucl. Phys. B38 (Proc. Suppl), 235 (1995). 

[24] L. Krauss, E. Gates and M. White, Phys. Lett. B299, 94 (1993); P. I. Krastev and S. 
T. Petkov, Phys. Lett. B299, 99 (1993); N. Hata and P. Langacker, Phys. Rev. D50, 
632 (1994). 

[25] S. P. Mikheyev and A. Yu. Smirnov, Yad. Fiz. 42, 1441 (1985); L. Wolfenstein, Phys. 
Rev. D17, 2369 (1985). 

[26] Frejus Collaboration, Ch. Berger et al., Phys. Lett. B227, 489 (1989). 

[27] NUSEX Collaboration, M. Aglietta et al, Europhys. Lett. 8, 611 (1989). 

[28] T. K. Gaisser, T. Stanev and G. Barr, Phys. Rev. D38, 85 (1988) and Phys. Rev. D39, 
3532 (1989). 

[29] G. L. Fogli, E. Lisi and D. Montanino, Phys. Rev. D49, 3626 (1994); CERN preprint 

25 



CERN-TH.7491/94. 
[30] L. V. Volkova, Sov. J. Nucl. Phys. 31, 784 (1980) . 

[31] M. Honda, K. Hidaka, and S. Midorikawa, Phys. Lett. B248, 193 (1990); H. Lee and 
Y. S. Koh, Nuovo Cim. 105B, 883 (1990). 

[32] E. L. Wright et al, Astrophys. J. 396, L13 (1992); M. Davis, F. J. Summers and D. 
Schlegel, Nature 359, 393 (1992); A. N. Taylor and R. Rowan-Robinson, ibid 359, 396 
(1992); J. A. Holtzman, and J. R. Primack, Astrophys. J. 405, 428 (1993). 

[33] J. R. Primack, J. Holtzman, A. Klypin, and D. O. Caldwell, UC Santa Cruz preprint 
SCIPP 94/28 and references therein. 

[34] J. Peltoniemi and J. W. F. Valle, Nucl. Phys. B406, 409 (1993) 

[35] J. Peltoniemi, D. Tommasini, and J. W. F. Valle, Phys. Lett. B298, 383 (1993); L. 
Bento and J. W. F. Valle, Phys. Lett. B264, 373 (1991). 

[36] D. O. Caldwell and R. N. Mohapatra, Phys. Rev. D48, 3259 (1993); 

[37] G. Raffelt and J. Silk, HEP-PH-9502306; G. M. Fuller, J. R. Primack and Y-Z. Qian, 
DOE-ER-40561-185. 

[38] D. O. Caldwell and R. N. Mohapatra, UC Santa Barbara preprint UCSB-HEP-95-1. 

[39] H. Minakata Tokyo Metropolitan University preprint TMUP-HEL-9502. 

[40] W. Pauli, on his famous letter presented to the Tubingen Meeting, in December, 1930 
introduces the neutrino as a "desperate remedy". 

[41] X. Shi, D. N. Schramm and B. D. Fields, Phys. Rev. D48, 2563 (1993). 

[42] V. Barger, K. Whisnant, and R. J. N. Phillips, Phys. Rev. D22, 1636 (1980); Phys. 
Rev. D43, 1110 (1991). 

[43] T. K. Kuo and J. Pantaleone, Phys. Rev. D35, 3432 (1987). 



26 



[44] A. De Rujula, M. Lusignoli, L. Maiani, S. T. Petkov and R. Petronzio, Nucl. Phys. 
B168, 54 (1980); H. Blumer and K. Kleinknecht, Phys. Lett. B161, 407 (1985). 

[45] E. W. Beier et al, Phys. Lett. B283, 446 (1992); K. Hidaka, M. Honda and S. Mi- 
dorikawa, Phys. Rev. Lett. 61, 1537 (1988); Phys. Rev. D44, 3379 (1991). 

[46] V. Barger and K. Whisnant, Phys. Lett. B209, 365 (1988) . 

[47] C. L. Smith, Phys. Rep. C3, 261 (1972). 

[48] G. L. Fogli and G. Nardulli, Nucl. Phys. B160, 116 (1979); M. Nakahata et al, J. Phys. 
Soc. Japan 55, 3786 (1986). 

[49] J. P. Revol, CERN-PPE/93-84 (1993), Presented at the "Rencontres de Physique de la 
Valle d'Aoste 1993", 8 March, 1993. 

[50] S. Mishra, private communication. 



27 



TABLES 

TABLE I. Characteristic values of L/E (m/MeV) for the different experiments considered in 
the analysis. 



Experiment 


(L/E) (m/MeV) 


Solar neutrino 


~ 10 10 


Atmospheric neutrino 


~ 10 2 


Bugey 


~ 20 


CDHSW 


~ 0.8 


E776 


~ 0.4 


E531 


~ 0.04 


CHORUS/NOMAD 


~ 0.02 


E803 


~ 0.01 


MINOS 


~ 70 
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TABLE II. Summary of the performance parameters for CHORUS and NOMAD. 





TV T S~\ TV T A T"~v 

NOMAD 


CHORUS 


Target 


Drift chamber panels 


Emulsion 


Mass (t) 


~ 3 


0.8 


Target thickness (in Xq) 


~ I 


~ 4 


Initial run period (yr) 


3 


3 


AT 1 £ T J. j_ • 

Number of CC Interactions 


f.f x 10 


5 x 10 


r identification 






Technique 


Pure kinematics 


Pure vertex 


Total pflfiripnrv in % fP/7f? y R fiM 


^ o 




Sensitivity: 






p 


1.3 x 1(T 2 


0.8 x 10~ 2 


P 

1 [IT 


2 x 10~ 4 


1.4 x 10~ 4 



TABLE III. Summary of the performance parameters for E803 and MINOS. 





E803 


MINOS 


Mass (t) 


~ 0.75 


10,000 


Target thickness (in Xq) 


~ 3 


Very Large 


Initial run period (yr) 


4 


2 


Number of CC interactions 


~ 9 x 10 6 


~ 10 5 


r identification 






Technique 


Pure vertex 


Disappearance experiment (ratio NC/CC) 


Total efficiency in % (EFF x BR) 


~ 10 


~ 70 


Sensitivity: 






p 


7 x 10~ 4 




P 


1.4 x 10~ 5 


0.012 
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TABLE IV. Summary of the performance parameters for NAUSICAA 





NAUSICAA (CERN) 


AT a TiriTr* A A /t — (TV T 1 T \ 

NAUSICAA (FNAL) 


Target 


Silicon/Graphite 


Silicon /Beryllium 


Mass (t) 


~ 1 


~ 2 


Target thickness (in Xq) 


~ 1 


~ 1 


Initial run period (yr) 


4 


4 


Number of CC interactions 


1.3 x 10 


2.4 x 10 


r identification 






Technique 


Vertex and kinematics 


Vertex and kinematics 


Total pffiripnrv in % fP/7f? y R fiM 




10 


Sensitivity: 






p 


7 x 10- 4 


7 x 10~ 5 


P 

1 flT 


2 x 10~ 5 


2 x 10~ 6 
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FIGURES 




FIG. 1. 90% CL exclusion contours for different reactor and accelerator experiments in the 
two-family mixing approximation: E531 — > v T apperarance experiment (solid line), Bugey 

"V^e disappearance experiment (dashed line), E776 — > z/ e apperarance experiment (dotted 
line), and CDHSW disappearance experiment (dot-dashed line). The allowed regions lie 

below the curves. 
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FIG. 2. (a) The 2a allowed regions in the Am\ T , sm 2 (29 flT ) plane for the atmospheric neutrino 
experiments for zero mixing er and allowed values of mixing e/i: Frejus (solid line), Kamiokande 
sub-GeV events (dashed line), 1MB (dotted line), and Kamiokande multi-GeV events (dot-dashed 
line). The arrows point towards the allowed regions for the different experiments. The thick solid 
line surrounds the 2a allowed region for the combination of all experiments. 

(b). Allowed region from the atmospheric neutrino deficit analysis together with the relevant 
constraints from the laboratory experiments for zero mixing er and allowed values of mixing e/x. 
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FIG. 3. (a) Same as Fig. 2. a for maximum allowed mixing er and allowed values of the mixing 

efi. 

(b) Same as Fig. 2.b for maximum allowed mixing er and allowed values of the mixing e/u,. 

(c) . Same as (b) but for maximum allowed mixing e/j, and corresponding allowed value for AMp M 
from the E776 experiment data. 
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FIG. 4. f e and spectra for the CERN beam (solid lines) and the Fermilab beam (dashed 
lines) as labeled. 
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FIG. 5. Accessible regions (90% CL) for the v e — > v T oscillation in the (AMp M , sin 2 (2# eT )) 
plane for the CHORUS/NOMAD experiments (fine solid line) and Fermilab experiment E803/the 
improved CERN experiment (dotted line). The dot-dashed line delimits the region accessible to 
NAUSICAA at Fermilab. Also shown in the figure are the region at present excluded by the Bugey 
data (shaded area) and the favoured range on AM^ M from dark matter considerations (solid 
horizontal lines). 
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FIG. 6. Accessible regions (90% CL) for the — > i/ T oscillation in the (AMp M , sin 2 (20 e(ti )) 
plane for an optimum value of the mixing er. The fine solid line corresponds to the CHO- 
RUS/NOMAD experiments and the dotted line corresponds to the Fermilab experiment E803 
and to the improved CERN experiment. The dot-dashed line corresponds to NAUSICAA at Fer- 
milab and the dashed line corresponds to MINOS. Also shown in the figure are the regions at 
present excluded by E776 data (shaded area) and the favoured range on AM|, M from dark matter 
considerations (solid horizontal lines). For comparison also the LSND data are shown. 
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FIG. 7. (a) Accessible regions (90% CL) for the — > z/ r oscillation in the (Am^ T , sin 2 (26> MT )) 
plane for zero value of the mixing er. The fine solid line corresponds to the CHORUS/NOMAD 
experiments and the dotted line corresponds to the Fermilab experiment E803 and to the improved 
CERN experiment. The dot-dashed line corresponds to NAUSICAA at Fermilab and the dashed 
line to MINOS. Also shown in the figure are the regions at present excluded by CDHSW data (dark 
shaded area) and the atmospheric neutrino analysis (light shaded area). 

(b) Same as (a) for maximum mixing er. Also shown in the figure are the region at present 
excluded by CDHSW data (dark shaded area) and the atmospheric neutrino analysis (light shaded 
area). These are the only constraints for zero or values of the e/i mixing sin 2 (20 e/1 ) <C 0.003. For 
maximum value sin 2 (2# e ^) ~ 0.003 the E776 limit (thick dot-dashed line) is relevant. 
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